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Combining the QCD dipole model with the Good and Walker picture, we formulate diffractive 
dissociation of a photon of virtuality Q 2 off a hadronic target, in the kinematical regime in which Q 
is close to the saturation scale and much smaller than the invariant mass of the diffracted system. 
We show how the obtained formula compares to the HERA data and discuss what can be learnt 
from such a phenomenology. In particular, we argue that diffractive observables in these kinematics 
provide useful pieces of information on the saturation regime of QCD. 
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\ I. INTRODUCTION 

o . 
o ■ 

. Hard diffraction has triggered a wide interest since its discovery in high energy virtual photon-proton scattering at 
HERA 0. It has been a major theoretical challenge to understand the observed large rate of diffractive events within 
. QCD, especially for large virtualities Q 2 of the photon (for a review, see Ref. Q). Among the attempts to describe hard 
f*^ ' diffraction, the dipole picture of QCD [jj 0] including unitarization corrections p| has proved particularly successful 
and natural in embedding both inclusive and diffractive observables in the same conceptual framework IE S ■ 

There are two distinct processes contributing to diffractive final states. A first one consists in the photon splitting 
in a qq pair, which scatters elastically off the target without any further radiation. The final state has typically a low 
invariant mass Mj, of the order of the virtuality of the photon. A second contribution is due to the qq pair interacting 
' through higher Fock state fluctuations, for example qqg, which go instead to a large-mass final state Mx ^S> Q- The 
£SJ , latter process is called diffractive dissociation of the photon. As there is a clear kinematical separation between the 
two processes, it makes sense to analyse these contributions separately. The present paper deals with the second one, 
04 , i.e. diffractive dissociation. 

To select this particular process one needs a large hierarchy between the mass of the diffractive final state and 
the virtuality of the photon. This means in practice that one is forced to relatively low virtuality scales Q 2 < 
1 GeV 2 , because of limited total center-of-mass energy available in the experiments. However, the process can still 
be computable in perturbative QCD provided that the energy is high enough to push the intrinsic k± of the partons 
£p-t inside the proton to large values k^^Q s ^> AqcTI • This feature is a consequence of the saturation of parton densities 
O j' and is now well-understood within QCD Q s is called the saturation scale, and is an increasing function 

(D , of the rapidity of the proton. 
i-C ■ The dipole picture is well-suited because color dipoles, which are two-body color neutral objets characterized by 
their longitudinal momentum, transverse size and impact parameter, are eigenstates of the QCD interaction at high 
energy. Both elastic scattering (and thus, from the optical theorem, the total cross section) and diffractive dissociation 
have a strai ghtf orward formulation in terms of dipoles. The latter is most simply obtained from the Good and Walker 
mechanism jlall3j . 

Our interest in diffractive dissociation is triggered by the wealth of new data that are being taken in the relevant 
kinematical r ang e at HERA. On the other hand, while some theoretical calculations are already available in the 
literature |l4Lll5j | . to our knowledge there has been no phenomenological analysis of this particular kinematical domain 
within QCD saturation models. Ref. 5] provides a model also for diffractive dissociation, but the approximations 
made there allowed only for high virtualities of the photon Q^>Q S , therefore the obtained formulae are unapplicable 
to the recent data at low Q 2 . 

In the following, we wish to clearly distinguish what is theoretically well under control and which are the model 
assumptions that have to be made to come to the description of the HERA data. In this spirit, we will start by 
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FIG. 1: One of the relevant graphs contributing to diffractive dissociation of the photon. The blob denotes the elastic double 
dipole-proton cross section. 

providing a complete derivation for the idealized process of diffractive dissociation of a photon on a large nucleus, for 
which the saturation corrections can be implemented in an easier way. This will be done in Sec. ^ Some analytical 
results are given in Sec. lIIII In Sec. lIVI we turn to phenomenology and explain how we can adapt the obtained results 
to photon-proton reactions at HERA. We compare our predictions to recently analyzed data. Our conclusions are 
drawn in Sec. E] 

II. DIFFRACTIVE DISSOCIATION OF A VIRTUAL PHOTON ON A LARGE NUCLEUS 

In this section, we compute the cross section for diffractive production of a hadronic system of large mass M\ 
in a highly energetic photon-nucleus collision. We limit ourselves to diffracted systems which have the following 
partonic content: qq or qqg. Higher Fock states are indeed not needed for large but still moderate values of Mj. The 
target is a large nucleus because, as will be discussed, complete and numerically manageable results can be obtained 
strictly speaking only in this idealized framework. We will also stick to the large- N c limit, and to the semi-classical 
approximation implied by the high energy kinematics. 

A. High energy kinematics and dipoles 

We go to a frame (the so-called dipole frame in which the photon and the nucleus have respective 4- momenta 
q 1 * = {q + 7 —Q 2 / (2g + ), X ) and — (0,p~, 0^). The total squared center-of-mass energy is given by s = 2q + p~ — Q 2 . 
Through an appropriate boost along the collision axis, the value of q + is chosen such that the whole diffracted partonic 
state \qq) or \qqg) is viewed as a quantum fluctuation of the photon which subsequently scatters off the nucleus. The 
nucleus contains all partonic fluctuations due to the remaining energy evolution as part of its wave function at the 
time of the interaction. This is possible in a frame in which the rapidity of the slowest final-state particle is close to 
zero. Of course, the observables will not depend on the frame, but our choice provides a physical picture which is 
particularly adapted to the study of unitarization effects. The event is diffractive if the nucleus does not break up, 
which requires the scattering to proceed through the exchange of a color-neutral object. 

It was proven [IcL ITU ] (for a review, see ^3) that in the dipole frame, the mean transverse momentum of the 
partons in the nucleus is shifted from Aqcd to a scale Q s (the so-called saturation scale) which increases with the 
rapidity of the nucleus, due to multiple scatterings and saturation of the parton densities [a 13 that occur to preserve 
the unitarity of the scattering matrix. If the energy is large, then the saturation scale can become high enough 
Q s 3> Aqcd to justify perturbative calculations for the diffracted partonic system, even if the transverse momenta at 
the level of the photon, initially of order Q, are not large. 



In the dipole frame, the time sequence of the interaction is the following: the photon couples to a qq pair, which 
then radiates (or not) a gluon before the whole system scatters off the nucleus through the exchange of a color-neutral 
object. The leading terms at high center-of-mass energy s come about when the ordering in time is strict, i.e. when 
the intermediate qq can be considered an almost on-shell onium state. The following factorization then holds: 

^£ = J d\ j\z(\r T (r,z;Q)\ 2 + K(^;Q)| 2 )^f , (1) 

where r is the transverse size of the qq pair, b its impact parameter, and z is the longitudinal momentum fraction of 
the antiquark. da^ff/d 2 b is the onium-nucleus diffractive cross section. The photon wave functions ifipi i\P L on a qq 
pair are given by [ijj 

\r T (r, z; Q)\ 2 = ^ £ e} {{z 2 + (1 - z) 2 )e 2 K 2 (e f \r\) + ™)K 2 {e f \r\)) 

3a ' (2) 
\r L (r, z; Q)| 2 = ^ £ e 2 4Q 2 z 2 (l - z) 2 ^ 2 ( £/ |r|) 

/ 

for a transversely and longitudinally polarized photon respectively. Here we sum over both polarizations, but in 
practice in the kinematical domain of interest (Q 2 not too large), only the transverse component will be relevant. In 
the previous formula, m/ is the mass of quark / and 



yj z(\ — z)Q 2 + to, 2 . (3) 



Thanks to the factorization Q), it will be enough to compute the cross section da^/d 2 b for diffractive dissociation 
of an onium. Using that formula, we will be able to relate it in fine in a straightforward way to the deep-inelastic 
scattering observables. 

Recalling that q + is the lightcone momentum of the incident photon, the 4-momenta of the partons in its wave 
function read 

h 2 h 2 h 2 
quark: (z q q + , r,fc g ) , antiquark: (z q q + ,- -,k q ) , gluon: (z s q , - + ,k g ) (4) 

ZZqq ZZqq ^Zgq 

where z q + Zq + z g = 1 and k q + kq + k g = from conservation of the 3-momentum. The transverse momenta k qi 
kg and k g are of the order of the external mass and/or virtuality scales which we assume all of the same magnitude 
Q 3> Aqce). We require furthermore the diffracted system to have a large mass in the sense Mx Q. Then the 
squared invariant mass of this partonic system reads 

M| = Q 2 (1 + 1 + 1) (5) 

\^q z q /- g J 

and gets large if either z q , Zq or z g are small. Due to the infrared singularity of QCD, z g is very small in a typical 
configuration, in which case M x — Q 2 / z g . This configuration with strong ordering of longitudinal momenta z g <t^. 
Zq, Zq is the dominant contribution in the kinematical regime of interest. Thus the mass of the diffractive final state 
is directly related to the longitudinal momentum fraction of the photon carried by the gluon. 

Let us now introduce some more kinematics. The difference between the rapidity of the onium and the rapidity of 
the slowest particle in the final state is denoted by yx ■ If (3 is the fraction of momentum of this particle with respect 
to the initial-state onium, then it is easy to see that yx — log(l//3). One also sees that (3 = z g for a \qqg) final state. 
The total rapidity difference between the photon and the nucleus is denoted by y = log(l/xBj), where xsj is the 
Bjorken variable. Then one introduces the rapidity gap variable 1 yp — y — yx and similarly to /3, the variable xp such 
that yp = log(l/xp). The relationship between these variables and the kinematics of the event reads 



Q 2 Q 2 x B j Q 2 + M 2 X 



Note that the saturation scale, which characterizes the transverse momenta at the level of the nucleus, depends on xr, and grows when 
Xf decreases. 



(a) (b) 



FIG. 2: (a) The qqg Fock state decomposition of the photon, (b) The same in the large- N c limit. 

The quantity which is measured and which will be computed is the differential cross section da 1 ^/ dMx or equivalcntly 
da^j d\og{\ / (3) . They are related by a simple kinematical factor: 

dg diff = 2M X da^ff 2 da^ 

dM x Q 2 + M| d!og(l//3) MJ»Q2 M x dlog(l//3) ' 1 ' 

The calculation of the cross section for diffractive production of a gluon in onium-nucleus scattering will proceed 
in two steps. In the next subsection we compute the elastic scattering amplitude for the onium on the nucleus up to 
order 0(a s ). Then we will deduce the diffractive cross section from the Good and Walker formula and with the help 
of the elastic amplitude. 

B. Elastic amplitude at next-to- leading order 

We only keep the leading terms in a 1/N C expansion. Then, as well-known, the partonic content of a color neutral 
object can be represented as a set of dipoles (see Fig. EJ. 

We denote by xi (a 2-dimensional vector) the transverse position of the dipole constituent i with respect to the 
center of the target, Xij = Xi — Xj, and Zi is the fraction of longitudinal momentum carried by constituent i. The 
normalization of a generic dipole state is given by the orthogonality condition 

(x, b, z\x', b', z') =8 2 {x- x')5 2 (b - b')S(z - z') (8) 

for two dipoles of respective sizes x, x' and center-of-mass coordinate vectors b, b' . In the following, unless it might 
hamper understanding, one does not keep track of the z-variables in order to simplify the notations. 

The physical dipole state |xoi, fr) p hys m the initial state of the experiment is dressed with all possible fluctuations. 
Its expansion in terms of bare Fock states reads 



Noi,&)phys = y/Z(x i)\x 01 ,b) + / d 2 X2dZ2^(x 2,Xi2,Z 2 )\xo2,Xi2,b) H (9) 



/phys 

where 4 > ( x o2, %12, ^2) is the probability amplitude that the dipole xqi split in two dipoles of respective sizes X02 and 
X12, which corresponds to a 1-gluon fluctuation of the initial onium. The gluon has longitudinal momentum fraction 
Z2- Note that \ip\ 2 is of order a s , and the higher Fock states not written are at least of order a 2 : this remark provides 
the systematics for this Fock-state expansion, seen as a perturbative expansion in the parameter a s . This expansion 
is relevant since a s is small because it runs with a scale between Q s and Q, which we have both assumed much larger 
than Aqci> Also, note that the 2-dipole state is not renormalized to this order, but would get renormalized from 
higher orders. 

The quantity that will appear in the calculation of observables is the squared wave function, summed over the 
polarizations and colors of the final-state gluon, whose expression can be found in 

<f>{x02,Xl2,Z 2 ) = ^ \ , 4>(X02,X 1 2 7 Z2)\ 2 = ~ 2^2 ' ( 10 ) 

polarization & color 

From the unitarity condition JSJ applied to the physical state |xoi, fr) p hys i n © ; one can compute the renormalization 
constant Z 

Z{x Q i) = l- I d 2 x 2 / dz 2 (f>{xo2,x 1 2,Z2) , (11) 
Jp 2 J/3 



where p cuts off the ultraviolet divergences and f3 the infrared ones. In practice, as it should be p disappears from 
physical quantities, and (3 is fixed by the kinematical limit. Note that Z represents physically the probability that 
the initial onium does not radiate any gluon. 

Let us denote by § the scattering matrix of the onium on the nucleus. Color dipoles form an orthogonal basis of 
the interaction eigenstates at high energy, thus the S-matrix element between two dipole states reads 

(x, b\S\x', b') = 5 2 (x - x')5 2 {b - b')S(x, b, x P ) . (12) 

The quantity \S(x, b, xp)\ 2 is the probability that a dipole of size a; at impact parameter b does not interact with the 
target. The xp-dependence of this quantity reflects the high energy evolution of the target. 
The S-matrix element between two double dipole states reads 

{x 2,xi2,b\S\x' 02 ,x' 12 ,b') = S 2 (x 02 - x' Q2 )S 2 (xi2 - x' 12 )S 2 (b - b')S(x m ,x 12 ,b,x r ) . (13) 

The notation is the following: S{xq 2i x\ 2 , b, xp) is the elastic diffusion matrix element for states formed of two adjacent 
dipoles of respective sizes xq 2 and xi 2 , and b = (xq + x\)j2. The quantity \S{xq 2 ,x\ 2 , b, xp)\ 2 is the probability that 
neither xq 2 nor x\ 2 interact, which under the assumption that the dipoles interact independently is the product of the 
probabilities that xq 2 and x\ 2 do not interact, \S(xq 2 , b — x\ 2 /2, xp)\ 2 \S{xi 2 , b — xq 2 /2, xp)\ 2 . The latter assumption 
is justified when the target is a large nucleus, since nucleons are color-neutral and independent. As furthermore S is 
known to be real at high energy, one can get rid of the squared modulus and one gets 

S(x 02 ,xi2,b, x P ) = S{x 02 ,b - xi 2 /2,x P )S{xi 2 ,b - x 02 /2,x ¥ ) . (14) 

Next one expresses the scattering matrix for the physical dipole \xqi, &) p hys as 

phys(xoi,b\S\x' Q1 ,b') phys = Z(x 01 ,b\S\x' 01 ,b') 

+ d 2 x 2 d 2 x' 2 dz 2 dz 2 %p*(x 2 1 x 12 ,z 2 )'ilj(x' 02 ,x' 12 ,z 2 ) (x 02 ,x 12 ,b\S\x' 02 ,x' 12 ,b') . (15) 

Using Eqs. I|11I12I14(I . and performing the integration over x' 01 and b' , one gets for the diffusion amplitude 

T(x x, b, (3, xp) = T(x i,b, x r ) 

a s N c f 1 dz 2 f 2 x^ 



2vr 2 



p z 2 



d x 2 2 — o - (^(^oi, b, x v ) - S(x 2,b-x 12 /2, x r )S(xi 2 , b-x 02 /2, x v )) , (16) 



where T stands for the matrix elements of T = 1 — § between bare states, and T are its matrix elements betweeen 
dressed states. The UV cutoff has been removed because when either xq 2 or x\ 2 tends to 0, then because of color 
transparency, the factor built from the S-matrix elements vanishes like x\ 2 (resp. x\ 2 ), and thus the integral over x 2 
is finite. The first term in Eq. H16[) is the Born term, the second one is the 0{a s ) correction to it. Physically, T is 
the amplitude for inclusive production of a forward gluon having longitudinal momentum larger than Pq + , when the 
total available rapidity is y = log(l//3a:p). 

At this point, it is interesting to note that the procedure described above can be iterated to get the Balitsky- 
Kovchegov (BK) equation ^]J. The latter gives the evolution of the dipole- nucleus S-matrix element for a dressed 
dipole S(xoi,b, f3,xp) when the total rapidity y — log(l/xBj) =log(l//3a;p) increases. Starting from the scattering of a 
bare onium at f3 = 1, Eq. H16J) gives the evolution of the amplitude T = 1 — S^xoi, b, /?, xp) as one decreases (3, i.e. as 
one gives a small boost to the onium while the rapidity of the target is kept fixed. The change of the amplitude comes 
from the fact that the onium can now be found with an extra gluon in its wave function: in the large N c limit, this 
means that the onium has a probability to be formed of two adjacent dipoles. Then under the assumption that the 
dipoles interact independently, for example with different nucleons of a large nucleus, this procedure can be iterated 
by simply replacing the S-matrix element for bare dipoles in the r.h.s. of Eq. Ijltjfl (for which (3 — 1) by the §-matrix 
element for the scattering of dressed dipoles at total rapidity log(l/2;2a;p) < log(l//3a;p), which reads 5(2:01,6, z 2 ,xp). 
One gets the following integral equation: 

S(x Q i , b, (3, xp) = S(x i , b, xp) 

+ 1 2 I / d 2 x 2 2 01 2 [S(x 02 , b- x 12 /2,z 2 , x P )S(x 12l b- x Q2 /2,z 2 , x r ) ~ S(x i, b,z 2 , x P )) . (17) 

^ Jf3 z 2 J X 02 X 12 V ' 



One now takes the derivative of both sides with respect to log(l//3). As we keep xp fixed, this is equivalent to taking 
the derivative with respect to the total rapidity log(l/x B j) and one gets 2 



dS(x 01 ,b,x Bj ) _ a s N c f 2 
d\og{l/x B] ) ~ 2tP " 



2 

x 2 2 01 2 (S(x 02 ,b- x 12 /2,x Bj )S(x 12 ,b- x Q2 /2,x Bj ) - S(x 01 ,b,x Bj )) , (18) 

X 02 X V2 



which is indeed the BK equation (see Ref. ^(| for a derivation closely related to ours). This equation resums the 
so-called "fan diagrams" . 

Coming back to Eq. I|16(l for T, the total cross section stems from the optical theorem, and the elastic cross section 
is the squared amplitude T: 



^ = 2T(x 01 ,b) and ^± = \T(x 01 ,b)\ 2 ■ 



(19) 



C. Diffraction from the Good and Walker picture 

An elegant and most straightforward formulation of diffractive dissociation is obtained by using the Good and 
Walker [l3] picture (see also [2(j and [jjj ) , according to which inelastic diffraction is proportional to the dispersion in 
cross sections for the diagonal channels of the scattering. The reason why we are relying on such a picture, developed 
in the context of early hadronic physics 0] and subsequently extended to the parton model pflj . is that the dipoles 
form precisely a complete set of eigenstates of the QCD interaction at high energy. The cross section for diffractive 
dissociation reads in this picture 



d?b 



f-= I d 2 x l m d 2 b'^ p ^ ys {x 01 ,b\T^\X){X\T\x' Q1 ,b') p ^ YS - / d 2 x' 01 d 2 b' phys {x i , b\T\x' 01 , &')p hy . 



(20) 



where X is any dipole final state, and the sum over X also contains an integration over phase space of the type 
J Y[ d 2 x a dz a . One denotes by (I) and (II) the (positive) terms appearing in this equation. 

The second term is the easiest to compute since it is simply da e \/d 2 b (see formula (2HJ)). Keeping only the terms 
relevant at order 0(a s ), one gets 



(") = = ^(soi, b, xp) + ^ I' — I d 2 x 2 ^k r {l - S(x Q1 ))(S(x 01 ) - S(x Q2 )S(x 12 )) 
d l b tH Jp z 2 J Xq 2 x1 2 



(21) 



Coming back to the first term in Eq. 112011 . the sum over X is restricted to single \x a b] b) and double \x ac , Xbc] b) dipole 
states, since we limit ourselves to the order 0{a s ) of perturbation theory. At high energy the interaction with the 
target cannot change the number of dipoles in the wave function, the first term thus gives 



(I)= / d 2 x' Q1 d 2 b'\Z(x 01 ) I d 2 x ab d 2 b"(x 01 ,b\^\x ab ,b")(x ab ,b"\T\x' 01 ,b') 

+ J d 2 x 2 d 2 x' 2 d 2 x ac d 2 x bc d 2 b" ^ (x 02 , xi 2 )^p(x' 02 , x' 12 ) 

x (xo 2 ,x 12 ,b\T^\xac,Xbc,b")(xac,Xbc,b"\T\x' 02 ,x / 12 ,b') \ . (22) 



Eq. fnji enables one to express Z(xqi) and Eqs. (|13I14|) to replace the matrix elements. Then one uses the (5-functions 
coming in Eqs. (|13I14J) which are due to the fact that the dipoles are eigenstates of the interaction, to perform the 
integrations over a;^, x' 2 , x ab , x ac , x bc , b', b" . One finds 



< l! = T2 (*oi) + ^ I' ^ \d 2 x 2 ^ r {i I S ( , lr2) SU l2} ) 2 - .1 Su<n)) 2 } 



(3 z 2 



,,2 

XnnXi 



(23) 



b 02 J '12 



2 We suppress the dependence upon the two variables /3, xp in Eq. 1181 because an inclusive observable is only sensitive to xgj = /3ip. 
Our derivation holds in a specific frame defined by the rapidity sharing between the onium and the target given by (3, xp, but of course 
the physical observable computed here is independent of this choice. 



Putting the two parts (I) and (II) together and taking the derivative with respect to log(l//3), one gets after some 
straightforward algebra 



da dm a s N c 



d 2 bd\og(l//3) 2ir 2 



d 2 x 2 2 01 2 (S(x 02 ,b- x 12 /2, x P )S(x 12 ,b - x 02 /2,x P ) - S(x 01 ,b, x P )) 2 . (24) 

X 02 X V2 

A few technical remarks are in order. Our result is identical to the one obtained in earlier derivations within different 
frameworks (semi-classical ^3, eikonal 0])- It is also consistent, in the weak interaction limit in which T = 1 — S <C 1, 
with previous calculations [2l| , and in particular with the computation of all Feynman diagrams contributing to the 
qqg final state [22j. 

Diagrams with gluons emitted in the final state are not included in our calculation since the radiative corrections 
we have concern exclusively the wave- function: no final state interaction can be accounted for in our calculation. This 
class of diagrams contributes to the diffractive final state, but they cancel out if one does not measure the transverse 
momentum of the gluon (see [lfiH^ V Note that we could not obtain more detailed properties of the final state within 
our derivation, but as demonstrated earlier, thanks to the selected kinematics we do not need them for our purpose. 

Using the factorization and the kinematical relation JTJl , one gets the cross section for diffractive dissociation of 
a virtual photon in the form 

da^ _ oJVc 1 / 2 / . - • . ,.,,,,:> 



dM x 7T 2 M x 



d 2 x i / dz (\ip^(z,x i;Q)\ 2 + \ipl(z : x i;Q)\ 2 ) 







..2 



! -oi 

:E 62 a; 12 



d X2 ~ji — — I d b(S(x 02 ,b- x 12 /2,x r )S{x 12 ,b- x a2 /2,x P ) - Sixoiib^xp)) . (25) 



Equation (|25|l has to be evaluated numerically for a general dipole-nucleus §-matrix, but there are two interesting 
limits that can be studied analytically, namely Q 2 3> Q 2 and Q 2 <C Q 2 - This is done in the next section. 

III. ANALYTICAL INSIGHT 

We will need a few general properties of the §-matrix element: 

S(r) < 1 (black disk limit) (26) 
\r\>l/Q s 

1 — S(r) ~ Ql 7 " 2 (color transparency) . (27) 

|r|<i/Q s 

Let us also recall that this theoretical study assumes enough center-of-mass energy so that Q 2 3> Aq CD in the dipole 
frame. This also means that the dipole sizes entering formula 1)24(1 . of order 1/Q S , are always smaller than the typical 
scale 1/Aqc;d of the variations of the §-matrix as a function of the impact parameter. Thus S(xo 2 ,b— x\ 2 ji) ~ 
S(xo 2 ,b) and S{x\ 2l b — xq 2 /2) ~ S{x\ 2l b). We will not have to take care anymore of the 6-dependence in the 
§-matrix elements, because it factorizes completely. 



A. Onium-nucleus cross section 



Small onium: collinear region 



We now analyse formula l|24|l in the limit x^Q 2 <C 1, £(3:01) — 1- The dipole splitting probability is singular for 



either |x 02 | < |x i| or |i 12 | < |x i| 
ordering relation implies \xi 2 \ ~ |xoi| 
can then be rewritten as 



Let us consider the first case, the second one being exactly symmetric. The 
The contribution of these dipole configurations to the cross section in Eq. J2U 



a s N c 
2ir 2 



d 2 x 2 — (S(x 02 )-l) 2 



(28) 



In this domain one has 1 — 5'(xo2) ~ Q 2x o 2 ^ so ^ e integral vanishes as Q s Xq^ . 

A second kinematical region which could bring large contributions to the integral over x 2 is |xo2| 3* |^oi|j which 
implies |xq 2 | ~ |xi2 1. Then the contribution of these configurations to formula (|24ll is 



a s N c 
2ir 2 



d 2 X2 ^f(S 2 (x 02 )-l) 



(29) 



(a) 



(b) 



FIG. 3: (a) The diffractive final state for Q 2> Q s - (b) The same in the large- N c limit. 



Using Eq. {2J|, 

one sees that this integral is of order Qlx^ and thus this kinematical region dominates over the 
previous one. It corresponds to configurations for which the qq pair is small and well separated from the gluon, or, 
in momentum space, the gluon has a small transverse momentum with respect to the one of the quark (see Fig. |SJ. 
This is the collinear limit. The lower boundary of the integral Q29J) can be put to zero since the region |xo2| < \xqi\ 
gives a subleading contribution. The result reads most generally 



The integral appearing here is proportional to Q 2 . The value of the associated proportionality constant has to be 
computed in specific models for S. 

The interpretation of the factor (1 — S 2 (x, b)) 2 is the following. S 2 (x,b) is the probability that a qq dipole of size 
x at impact parameter b does not interact with the nucleus. The elastic cross section reads in this case do~ e \/d 2 b = 
(1 — S(x, b)) 2 . In the same way, S (x, b) would be the probability that two independent but superimposed dipoles do 
not interact, and da e \/d 2 b = (1 — S 2 (x, b)) 2 would be the corresponding cross section. So the quantity (1 — S 2 (x, b)) 2 
can be understood as the elastic cross section for an octet dipole, because in the large N c limit the latter is represented 
by two independent lines of indices. 



We now turn to the limit x^Q 2 3> 1. The discussion is more subtle in this case because unlike the previous case 
one cannot identify a typical dipole configuration which gives the dominant contribution to the cross section. Let us 
split the integration domain into several relevant subdomains as shown in Fig. 0] First consider the disks T>q and T>\ 
of radius 1/Q S around xq and X\ respectively, inside of which S(x\-i) ~ S(xoi) and S(xo 2 ) ~ 1 (resp. S{xq<x) ~ S{xq\) 
and S(x\i) ~ 1). Next define the domain X for which either \xqi\ or |xi2 1 is smaller than |xoi| (the intersection with 
T>q and T>\ is excluded). The remaining part of the 2-dimcnsional plane is denoted by £ . 

In the latter integration region £, necessarily S(xq 2 )S{xx2) <C S{xq\) but in I, the ordering relation depends on 
the model for S. Taking the relevant approximations in these different domains, the integral in Eq. (|24|l then reads 



The first term goes like S 2 {xq\), and the second term is enhanced by a factor of \og{Q 2 x 2 1 ). As for the third term, 
no completely general statement can be made. If S(xq2)S(x\2) <C S(xqi) everywhere in domain X, the third term 
contributes as much as the second term. Then the cross section goes like S 2 (xm ) logfccm Q 2 )- This is the case for 
color glass condensate for which S(xqi) ~ exp(— A log 2 (xQ 1 Q^)), with A ~ 0.2 j2J,|2||. If S^xc^S^x^) > -S'(xoi) 
instead in some region inside X (and in particular at the point |xq2| = |xx2| = |xoi|/2), the third term contributes 
as S 4 (xoi/2)/xq 1 , and can even dominate over the two other ones. This is what happens when S is given by the 
Golec-Biernat and Wusthoff model. Thus we see that the diffractive cross section behaves qualitatively differently for 
the various models for the dipole-nucleus S-matrix. To summarize, in these two kinds of model one finds the following 
estimates for the asymptotic behaviors: 




(30) 



2. Large onium: saturation region 



(xoi) / d 2 x 2 f^-(x 02 )V + S 2 (x 01 ) / d 2 x 2 ^V+ / d 2 x 2 ^ 

JT>a+T>i \CX02 / j£ X 02 X 12 Jl X Q2 X \ 



(S(xo2)S(x 12 ) - S(xoi)) 2 • (31) 




for CGC, 



for GBW, up to a 0(1) factor. 



(32) 



FIG. 4: The transverse plane and the different integration domains T>o, T>\, T and £ . The crosses denote the positions of the 
quark and the antiquark, xq and x\ resp. 



For GBW-like models, 0(1) factors are cumbersome to get. 

The onium cross section is anyway strongly damped in this region of large x^Q 2 because S(xqi) goes to in this 
limit, much faster than 1 / (\og(xg 1 Q 2 )) in all available models. This is the black disc limit, in which the elastic cross 
section is half the total one, and thus the diffractive cross section is zero from the inequality 



which directly follows from unitarity |2(| . 

To summarize, we see that the diffractive cross section increases with |xoi| like Xq X at small x^Q 2 and decreases 
strongly at large x^Q 2 , so there must be a maximum for x\ x ~ For a realistic target, the saturation scale is a 

decreasing function of the impact parameter. For a well-chosen value of xoi, it can happen that Q 2 s (j))x\ x 3> 1 for small 
b, and Qg(b)xgi -C 1 for large b (this is always the case). In such a situation, the main contribution to the diffractive 
cross section for an incident dipole of size \xqi\ comes from a zone in impact parameter which has the shape of a 
ring of radius bo such that Q 2 (bo) ~ l/^oi- We have recovered the well-known fact [2(J that diffractive dissociation 
is peripheral at variance with elastic or total scattering. However, the "periphery" of the nucleus is denned with the 
help of the saturation scale and not with Aq^tj at variance with soft processes. 

The features just found are illustrated on Fig. for the Golec-Biernat and Wiisthoff model, and for the color glass 
condensate in the large Xq^Q 2 limit. 



Before we turn to the comparison to the data, we wish to study some analytical properties of the cross section Q25JI 
in the two different limits Q 2 ^> Q 2 and Q 2 <C Q 2 - We focus here on transversely polarized photons. To study 
asymptotic behaviors, it is enough to replace the photon wave function that appears in Eq. (|25[) by the following 
approximation: 



da dig < 1 dvtat _ da e \ 
d 2 b ~ 2 d 2 b ~ ~d?b 



(33) 



B. Limiting behaviors of the photon- nucleus cross section 



\i/j^(x 01 ,z;Q)\ 2 




]T e){z 2 + (1 - zf)-^Q{l z(l z)Q 2 x 2 01 ) . 



01 



(34) 



10 




x 01 Q s 

FIG. 5: The diffractive cross section compared to the approximations in extreme kinematical regimes. 



3a fi 



(35) 



Let us first treat the integration over z in Eq. 1)25(1. There are two cases to distinguish. If Q 2 x\ x < 4 then the 
O-function in (|34(l can be ignored and 

/ dz |^J(xoi, z; Q)\ 2 = -^j-y2 e2 f-r- ■ 

Jo n f x oi 

If Q 2 x 2 ll >4 instead, then the leading contribution, obtained for Q 2 x 2 )1 ^>4, comes from the endpoints 

/•l/(Q 2 *oi) /•! 

dz\ip^,(x i,z;Q)\ 2 = / dz|-0y(a;oi,z;Q)| 2 + / |^(x i, 2; Q)| 2 



(36) 



This is the well-known aligned-jet configuration |26| . 

As for the case Q 2 3> Q 2 , the integral over xoi can be split in 3 integration regions [0,4/Q 2 ], [4/Q 2 , l/Q 2 } and 
[l/Q 2 ,oo]. Taking the relevant approximations for the wave function 1(35(1 . ((36(1 and for the dipole 8- matrix ele- 
ment (I3UII . ((32(1 in each of these domains, one sees easily that region [A/Q 2 ,l/Q 2 ] dominates. Eq. I(25|) reduces to 



dal m 6a em N c 1 



/ ^log|/$(l-^(-,6,^)) 2 ■ (37) 



This limit is consistent with previously derived results, as shown in the appendix. 

As for the case Q 2 <§C Q 2 , one now splits the integral over xqi in [0, l/Q 2 ], [l/Q 2 , 4/Q 2 ] and [4/Q 2 ,oo]. This time, 
the first domain gives the main contribution which reads 

da diS a em a s N cS -^ 2 1 f 2 f d 2 x 2 2 ,o Q \ 



14 



12 



10 



da diff /dlog(1/p)/xg r 
-d/dlog Q 2 (xoJgdz |y 7 | 2 ) for Q 2 =10 1 GeV 2 
} Q 6z |v|/ y | 2 forQ 2 =10 1 GeV 2 ■ 
Q 2 =10 2 GeV 2 
Q 2 =10 4 GeV 2 ■ 




FIG. 6: The dipole cross section and the photon wave function for three different values of Q 2 . The former is divided by the 
squared size of the dipole, and the latter is multiplied by the same, so that the convolution of these two quantities is equal 
to the convolution of the dipole cross section and photon wave function. The saturation scale Q 2 is set to 10 GeV 2 . For 
Q 2 — 10 GeV 2 , the derivative of the wave function with respect to log Q 2 is also shown. 



and which is independent of Q 2 and Q 2 (recall that the integral over the vector x is proportional to Q 2 S ). This can be 
illustrated graphically. As one goes to smaller Q 2 , Xq X J dzi\ip 1 (xoi, z\ ;Q)\ 2 stays constant for lower values of Xq X , 
and develops a plateau (2^, as seen on Eq. and on Fig. This plateau contains configurations for which the 
longitudinal momentum of the photon is equally shared between the quark and the antiquark. On the other hand, 
ticr ( jiff/dlog(l//3)/xQ 1 is roughly a constant for x\ x Q 2 <C 1 and drops rapidly to zero for x\ x Q 2 s ^> 1. Hence once 
Q 2 is small enough to allow the plateau of the photon wave function to extend up to dipole sizes of the order of 
1^01 1 ~ 1/Qs, the cross section does not increase anymore when Q 2 decreases, and tends to a constant. The precise 
value of the constant is model-dependent. 

These features are specific to the transverse photon and to the purely diffractive cross section. 



IV. PHENOMENOLOGY 



We are going to apply the formalism developed above to deep-inelastic scattering at HERA. Our starting point is 
Eq. 



A. From a nucleus to a proton target 



To go from scattering on a nucleus to scattering on a proton, one has to take some model assumptions for the 
non-perturbative physics inherent to the target and that are not under theoretical control. 

First, we assume that the incident dipole is always much smaller than the target proton, which is reasonable for 
Q 2 not too small compared to Aq CD . We also consider a cylindrical target, i.e. we neglect all dependence upon 
impact parameter. As far as we do not consider observables which would depend on the momentum transfer to the 
proton, this assumption only hampers predictivity for the global normalization, and brings a significant technical 



simplification. 

Although we are now dealing with a proton, we still keep the hypothesis of independence of dipole interactions in 
the derivation of Eq. I)25[l. which was a priori justified for a large nucleus only. This could be a problem for Q 2 < Q 2 . 
But it has recently been shown (see [2^, 113 and |2j| for the theoretical justification) that the proton can be seen 
effectively as formed of color neutral domains of typical size 1/Q S . Thus deep in the saturation regime Q 2 <C Q 2 , the 
proton is not different from the nucleus, in the sense that the dipoles in the projectile interact at most once with each 
of these color independent domains in the proton. 



B. The saturation model 

In practice, we choose to stick to the Golec-Biernat and Wiisthoff (GBW) model @ for the dipole-proton S-matrix, 
which reads 

S(r, b, x P ) = 0(R P - b)e- Q2 ^ r2 l A + 6(6 - R p ) , (39) 

where Q 2 (xp) — (x/xo)~ x in units of 1 GeV 2 and R p is the radius of the proton, related to a normalization parameter 
in GBW model by R 2 = ao/(2n). xo, A and a were fitted to the data for the total structure function F 2 in Ref. 
we just take over the found parameters Co = 23 mb, A = 0.288 and xo — 3.04 x 10~ 4 . The integration over the impact 
parameter that appears in Eq. I|25|l can be performed, and yields a factor ao/2. In addition, in the GBW model, 
three light quarks were considered, and a mass of 140 MeV was assigned to them, to ensure a sensible extrapolation 
to photoproduction. We take over this feature to our model. 

C. A comparison to the data 

We find useful for numerical evaluation to rewrite the onium diffractive dissociation cross section as 
da diS 2a s C F r ^ f 1 ^ 1 1 



d 2 bd\og{l/(3) TT 2 7 J l-A! + 2AiA 2 ^/\ 2 (1- A 2 )(l + AiA 2 )(l-Ai + AiA 2 ) 

x |i-E 2 (a; i;Ai;l-A 1 + 2A 1 A 2 ) + A 1 E 2 ( 2 ; 01 ;l/Ai;l/A 1 -l + 2A 2 )| , (40) 

where T,(xoi] a;b) = S(axoi)S(bxoi) — S(xqi). Through an appropriate change of variable, we have mapped the 
complex plane into the finite domain [0, 1] x [0, 1]. Furthermore, all the factors of the integrant are not more singular 
than l/\/A for A — > 0, which is integrable. So this formula does not require numerical cancellations between large 
terms, which would result in large errors. Note also that the obtained formula is quite simple, and this feature might 
be related to the conformal invariance of the dipole splitting kernel. 

Putting everything together, the formula that we have to evaluate numerically is 

d^s^^sNc^o r dxQlXQ1 f dz(W T (x i, z; Q)\ 2 + |^(x 01 , z; Q)\ 2 ) 
dM x 7T M x Jo Jo 

f 1 1 

x / d\idX 2 



[o,i]x[o,i] l-Ai + 2AiA 2 ^A 2 (l-A 2 )(l + AiA 2 )(l-Ai + AiA 2 ) 

x |^-S 2 (s 01 ;A 1 ;l-A 1 + 2A 1 A 2 ) + A 1 E 2 ( a; oi;l/Ai;l/A 1 -l + 2A 2 )| . (41) 

where £ is constructed from the §-matrix element (13911 . and the photon wave functions are given by Eq. with 
flavors and quark masses chosen as explained in Sec. IIV Bl 

The data do not distinguish between elastic and diffractive contributions, so we have to add a component in which 
the final state is a qq pair. This component is not of direct relevance for us as it gives a non-ne gligi ble contribution 
for large (3 ~ 1 only. As it has been extensively studied in the literature within different models EHL , we just 

reproduce the parametrization obtained in Ref. [B| to estimate the importance of this contribution. 

A priori, there is no free parameter. However, by taking a s at its value at a phenomenologically realistic scale 
Q 2 ~ Ql ~ 1 GeV 2 , we overshoot the data by a factor 1.5 — 2. We must set a s to a lower value a s = 0.15. 

However, the fact that we do not predict correctly the global normalization does not come as a surprise. Indeed, 
first, the (non-perturbative) assumption of a cylindrical target is certainly not realistic. Second, we took the size of the 
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FIG. 7: Comparison of the model with the data. The dashed line is our computation of the inelastic component i|4 ipi . The 
dotted line is the elastic component, which is taken from Ref. The full line is the sum of these two components, which for 
high diffracted masses, reduces to the inelastic component. The data points are taken from [33| . 



cylinder over from the normalization of the cross section fitted to total cross sections. Because diffractive dissociation 
is resonant with the saturation scale, the impact parameter region which contributes to the cross section should be 
effectively smaller 3 , which would go in the right way as far as the normalization is concerned. 



3 We thank Edmond Iancu for having drawn our attention to this point. 



Our predictions are presented on Fig. [7] As our model is established in the limit of small values of j3, we choose 
to restrict our calculation to (5 < 0.2. This explains why our curves do not extend to very large Q 2 , especially at 
Mx = 5 GeV. We see that we get a good agreement with the data. 

D. Future prospects on phenomenology 

On Fig. we see that the cross section tends to a constant at low Q 2 , as anticipated in Sec lIII Bl This was shown 
to be due to the fact that Xq X J q dz \?Pt\ 2 a l so tends to a constant at low Q 2 , see Fig. which is unavoidable when 
the initial state is a transversely polarized photon. The consequence is that the high sensitivity of these diffractive 
observables to the exact form of the §-matrix in the saturation region is somewhat spoiled by this smearing. 

However, this property is not true for a longitudinally polarized photon. Alternatively the plateau at small x^Q 2 can 
also be eliminated by considering the derivative of the cross section with respect to logQ 2 , namely daLg/dMxd log Q 2 
|27|. see Fig. EI A measurement of this quantity would help the understanding of the S- matrix in the saturation regime, 
by enabling a direct scan of this interesting region. A more detailed study will be provided elsewhere. 

V. CONCLUSION 

Diffraction is a good place to study saturation effects in deep-inelastic scattering. It had already been shown that 
quasi-elastic diffractive scattering, like vector meson production, enables to measure how far one is from the unitarity 
limit and gives a handle on the impact-parameter dependence of the saturation scale |34ll35j]. In this paper, we have 
provided a theoretical tool to study "true" (inelastic) diffraction in the saturation regime. We have formulated the 
diffractive cross section for a general model for the dipole-proton S-matrix, and we have rederived an already known 
formula in an elegant way using the Good and Walker picture. When we take for S the Golec-Biernat and Wusthoff 
model, we obtain parameter- free predictions (up to the global normalization) for the observable da/dMx currently 
measured with good precision at HERA. An important highlight of our analysis is that diffractive observables can 
help to discriminate in a unique way between the predictions of different models in the saturation region. This point 
deserves more studies, that we leave for the future. 

On the theoretical side, the formalism used here based on the dipole model could be easily generalized to an 
arbitrary number of gluons in the final state. This program has already been explored in Refs. Q, HfJ but no 
unitarity corrections were taken into account. Our method allows to take into account these corrections. Whether a 
simple evolution equation in f3 would be found is not clear. Anyway, the available energies at present colliders do not 
require yet such higher Fock states, nor the complete resummation of all leading logs, i.e. of the powers of a s log(l//3). 

On the phenomenological side, it would be worth to take into account more accurately the dependence upon impact 
parameter, which is of great importance in the whole discussion of saturation Therefore, one could replace the 

GBW ansatz for the dipole S-matrix by Kowalski-Teaney recent model, which takes into account the transverse 
profile of the proton. This replacement might be necessary in order to describe accurately enough the high precision 
preliminary data, and we would have a prediction also for the normalization of the cross section instead of tuning it 
as we have done here. Furthermore, more exclusive observables like da/dMxdt, which will be measured in the future 
could be predicted. 

APPENDIX A: HIGH MASS DIFFRACTION FROM THE COLLINEAR LIMIT 

In this appendix, we compare our approach (see Eq. I|37|)) to the impact factor used by Golec-Biernat and Wusthoff 
for qqg final states. The latter is formally valid only for Q 2 ^> Q 2 , i.e. in the collinear limit. The two calculations 
should however lead to the same result when both Q 2 3> Q 2 and (3 -C 1, which is the so-called double logarithmic 
limit. 



According to Ref . |a, |32] , the diffractive structure function reads 

da^ _ 81a em /3Mx \- 2 a» f 1 dz ( f f3\ 2 ( f3\ 2 \ z 
dM x 6A^BQ 2 (Q^+M%)^ 6f 27r J p z [\ z) + \z ) j (1 - z) 3 



d 2 k 



A 4 log ( (1 ^ Q2 \ 6((1 - z)Q 2 k 2 ) J d 2 rd 2 r'e^ r -^a gg (r,x P )a gg (r\x P ) 



x E(^ mn - - 2 ^)^ Ur- Z kr ) K2 Uih kr ' ) (A1) 



where a gg is the adjoint dipole-proton cross section, and B is the exponential t-slope of the differential cross section: 

e - Bt . (A2) 



d(J dm _ d<J diS 



dt dt 



One can extract the leading contribution to this cross section in the /3 — ► limit, ft comes from the integration 
region z ~ /3. Using i^C^) ~ 2/a; 2 and performing the sum over helicities m,n, one gets 



dal m _ 81a em /3 ^ 2 a s f 1 dz ( ( [3\ 2 //3\ 2 \ f d 2 k fQ 2 \^ t/) 2 ,2 



dM x 8it 3 BQ 2 M x ^ 1 2tt J z 2 \ \ z) \z J ) J (2tt) 2 b \k 2 



d 2 rd 2 r' ' <■•■ -'V-' 



^ (r - r 'V gg (r,x P )a ffg (r' 7a ; P ) (2^^- - (A3) 



The integral over z is now factorized and gives just a factor 2/3/3 in the limit j3 — ► 0. The integral over A; can be done 
analytically 

log (|!) e *-) 6 «3= - *■> - i '-^r'" . (A4, 

One performs the change of variable R — Q s r, X = Q s (r — r'). The following expression results: 

dal s _ 27a em ^ 2 a s f „. 1 - J ((Q/Q.)I*I) 



where we have defined 

f{\X\) = Q 2 s fj d<P J R2( x* x) 2 *M x P )a gg (R-X, x r ) ^ { %^Zxy - 1 ) ' (A6) 

and 4> is the angle of vector X, a is the dipole cross section after change of variable. /(|A'|) is a smooth regular 
function of \X\ having finite value at X = and decreasing like 1/|A| 2 . Let us choose a number A > 1 such that 
I Jo(A)| <C 1. The integral over \X\ can be decomposed as follows: 



/(0) 



XQ °' Q d \x\ l — M(Q/Q s )\x\)f(\x\) , /- 1 , f°°d\x\f(\x\) 



x\ /(o) 

JXQs/Q \x\ /(o) A |x| /(o) 



(A7) 



As the integrand is bounded by a number proportional to Q/Q s , the first term in the square brackets tends to a constant 
for Q 3> Q s . The third term is also a number since /(|X|) decreases sufficiently quickly to ensure convergence of the 
integral. The leading contribution to the second term is log(Q/Q s ). The coefficient of this logarithm is the value of 
/ at zero, namely 

/(0) = 2. / d 2 r^^l . (A8) 



Putting everything together back into Eq. I|A5() . one finally gets 



1 v-^ 2 a s , Q 2 f 



dal m = 27a em f Vji2 a s ^Q 2 f ^^a 2 g (r,x P ) 

dM x ~ Q 2 



d 2 r " v ; r ' . (A9) 



In order to compare with our approach, one has to recall (see Sec II V A|) that in the large- N c limit, a gg is given by 

a gg (r,x r )=2 [ d 2 b(l~S 2 (r,b,x P )) = a Q (l - S 2 (r,x P )) , (A10) 



where S is the 3 — 3 dipole S- matrix element. Replacing Eq. I|A10|) into Eq. (|A9|) . the ratio between (|A9|I and 1371) 

reduces to a constant: 

. (All) 

The fact that the results match only up to a factor is to be traced to the different treatments of the integration over 
impact parameter in the two cases. 
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